Magnetic groundstates in a correlated two orbital Hubbard model 
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We examine the orbital and magnetic order of the two orbital Hubbard model within dynamical 
mean field theory. The model describes the low energy physics of a partially filled e 9 -band as can be 
found in some transition metal compounds. The model shows antiferromagnetic as well as ferromag- 
netic phases. For stabilizing ferromagnetism we find that Hund's coupling is particularly important. 
Quarter filling represents a very special situation in the phase diagram, where the coupling of spin, 
charge, and orbital degrees of freedom are involved. Exactly at quarter filling we find a metal in- 
sulator transition (MIT) between two almost fully polarized ferromagnetic states. This MIT can 
be tuned by changing the local interaction strength and seems to be a first order transition at zero 
temperature. Apart from these ferromagnetic states we were also able to stabilize antiferromagnetic 
and charge ordered phases at quarter filling, depending on the interaction parameters. 

PACS numbers: 71.10.Fd,71.30.+h 
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I. INTRODUCTION 

Magnetism is still a very important topic in condensed 
matter physics. Analyzing the elements involved in mag- 
netic materials, one has to conclude that the existence 
of magnetism is intimately connected to the presence of 
partially filled d- or f-shells. It is thus not surprising 
that understanding the influence of these shells on the 
low temperature physics is crucial for a proper descrip- 
tion of magnetic materials. 

Since d- or f-shells are typically more strongly localized 
than the s- or p-shells of simple metals or semiconduc- 
tors, they often are subject to strong correlation effects 
making their theoretical study a challenging task. In 
this context very interesting materials are the transition 
metal oxides^. One prevalent lattice structure of tran- 
sition metal oxides is the cubic pcrovskite, in which the 
transition metal atom sits in the center of an oxygen octa- 
hedron. The states relevant at low temperatures are the 
d-orbitals of the transition metal atom. Due to the cubic 
crystal symmetry the d-orbitals split into a threefold de- 
generate tig- and a twofold degenerate e g -band, which for 
the coordination present in the pcrovskite structure has 
the higher energy compared to the i2g-band. Thus not 
only strong correlations but also orbital degeneracy plays 
an important role in the physics of these compounds. 

To give a specific example for a transition metal ox- 



or more pre- 



ide, let us consider the Manganites^i™^, 
cisely Lai-^Ca^MnOa. Manganites became famous for 
their colossal magneto resistance (CMR). Besides this 
particular feature they show a very rich phase diagram 
with different magnetic and orbitally ordered phases. In 
Lai_2;Ca x Mn03 one has to distribute 4 — x electrons 
per site to the d-states of the manganese according to 
Hund's rules. Thus, the electronic configuration in this 
compound can be modeled by a partially filled e s -band 
close to quarter filling and a half filled tig-hand which 
couples via Hund's coupling ferromagnetically to the e g - 
electrons. The hopping between the ^-states is very 



small and thus the ^-states are often modeled as lo- 
calized 5 = 3/2 spins. Besides this electronic part the 
lattice degrees of freedom and especially Jahn- Teller dis- 
tortions are important to correctly describe the physics 
of Manganites^. 



The ultimate goal surely is to theoretically describe the 
properties of materials like Lai_ x Ca x Mn03 including all 
the degrees of freedom mentioned before. However, at 
least equally important for a proper microscopic under- 
standing of the physics is to disentangle the contributions 
of the different degrees of freedom and identify the indi- 
vidual influence, in particular to what extent a certain 
degree of freedom is responsible for effects or just follows 
the leacU^. For this reason we now leave the special topic 
of Manganites and focus on the role of the electronic de- 
grees of freedom, in particular the role of interactions 
on the magnetic and orbital properties of a degenerate 
e s -band. We thus ignore effects of the i2 g -spin and the 
lattice in the following. Based on such an investigation, 
one can later include further features, like different band 
widths of the eg-ban d 11 : 12 ' 13 : 14 ! 15 , or additional degrees of 
freedom, for example the ^g-spin (or band) respectively 
the strong coupling to the lattice, step by step, thereby 
properly identifying for which particular features they are 
actually responsible. As we will discuss, already the sim- 
plified situation with only an e 9 -band present shows very 
complex ground state properties, involving the coupling 
of charge, spin and orbital degrees of freedom. 



The article is organized as follows: After this introduc- 
tion we establish the model and discuss the methods we 
have used for solving it. In the results section we first 
analyze the situation for a band filling between quarter 
filling and half filling. In the second part of the results 
section we particularly address quarter filling represent- 
ing a very special point in the phase diagram, at which 
different ordered phases compete with each other. 
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II. MODEL AND METHOD 



We model the e g -band as a two orbital Hubbard 



model^OLisds^ 

2 2 

ij,(T 771 — 1 2,cT m — 1 



- 2 J Si t iSi t 2 



m=l 



+(E/' - J/2) y^(nj,-|-,i + n ii j, i i)(n ii -|' i 2 + "1,1,2) 



Here i, j label the lattice sites; to = 1, 2 is the orbital and 
(7 the spin index. Thus, cj ^ m creates an electron at site 
i in orbital m with spin a. As usual, [i represents the 
chemical potential with n = c^c being the density opera- 
tor. Note that, in the spirit of the philosophy discussed 
in the introduction, we assume the same band structure 
and -width for both orbitals in the e s -band and also do 
not include orbitally off-diagonal hopping. 

The two particle interaction is parametrized as in 
Ref. [10, but we neglect the two particle hopping term 

J rc !,T ,m c !a,m c U,™ c i,T ,rr^- Tnis is done mainly for numer- 
ical reasons, as otherwise one cannot introduce a con- 
served orbital quantum number. This of course intro- 
duces an additional approximation; however, previous 
studies indicated that this particular term is of minor 
importance, at least for ferromagnetic Hund's exchange, 
while the inclusion of the rotationally invariant spin ex- 
change appears to be crucial 2 -^. Therefore, in our model 
the two particle interaction is given as a local intra-band 
interaction with amplitude U , a local inter-band inter- 
action with amplitude U' — J/2 and a Hund's coupling 
with amplitude 2J between the spins Si and S2 of the 
two orbitals. 

In order to investigate the possible magnetically, 
orbitally and charge-ordered phases of the two or- 
bital model we use the dynamical mean field theory 
(DMFT) 22 ' 23 . As this mean field theory properly in- 
cludes the local dynamics due to electronic correlations, 
it accounts for such subtle effects like crossover between 
itinerant and localized order— and strong reduction of 
transition temperatures. As it is well-known, the lattice 
structure enters the DMFT calculation only via the non- 
interacting density of states (DOS). Since we are inter- 
ested in general qualitative aspects of ordering phenom- 
ena in the two band model and not in the description 
of a particular material, we can use this property and 
choose a numerically convenient form of the DOS. As 
has been discussed extensively^ 2 -, the semi-elliptic DOS 
obtained from a Bethe lattice with infinite coordination 
number is indeed a numerically convenient choice, which 
also allows to address the ordered structures we are in- 
terested in, namely antiferromagnctic Neel, homogeneous 
ferromagnetic, as well as charge and orbital order. More 
complex structures, like for example ferromagnetic order 



in one direction and antifcrromagnetic in another— are 
excluded deliberately. 

The DMFT self-consistency maps the lattice onto an 
effective impurity Anderson model, which in the present 
case becomes a two orbital Anderson model. We solve 
this Anderson model with the Numerical Renormaliza- 
tion Group (NRG) 26 ' 27 , which allows us to calculate 
properties for wide parameter regions at T = 0. To re- 
liably calculate spectral function with NRG we use the 
complete Fock space metho d 28 ' 29 . It must be noted that 
the two orbital model is an extreme case for calculat- 
ing spectral functions within the NRG. The calculation 
of one spectral function requires approximately 15GB 
shared memory and several CPU-hours. As discretiza- 
tion parameter within the NRG we used A = 2, keeping 
up to 5000 states per NRG step in a typical calculation. 



III. RESULTS 

A. Filling range 1 < (n) < 2 

Let us begin the discussion of our numerical results by 
an overview of the magnetic phase diagram as function 
of filling for 1 < (n) < 2. The calculations were carried 
out using a fixed local intra-orbital interaction U/W = 4. 
Here, W is the bandwidth of the semi-elliptic DOS, which 
will be used as energy scale throughout this article. This 
strength of the interaction is a good guess for transition 
metal oxides. As in the one orbital Hubbard mode l 30 ' 31 , 
we also observe that for strong local interaction and half 
filling the physics is dominated by an antifcrromagnetic 



0.9 
o.s 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 




B 




0.8 1 1.2 

Occupation <n> 



1.6 



FIG. 1: (color online) Ferromagnetic polarization for J = W 
as function of filling at U = AW, V = U - 2 J and T = 0. For 
fillings n > 1.5 we observe signs of an incommensurate spin 
density wave extending to half filling. The phase diagram was 
created by fitting a smooth surface through approximately 
40 inhomogeneously distributed data points. Therefore the 
location of the phase boundaries are only meant as rough 
sketches. 
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insulator originating from super exchange 3 ^. The Neel 
temperature does not depend on J or U' within the tem- 
perature resolution given by the NRG, as long as U is the 
dominating interaction. The antiferromagnetic phase can 
be doped resulting in an incommensurate spin density 
wave away from half filling, which can extend to occupa- 
tions (n) ?s 1.5. This is again the behavior expected for 
strong local interaction 3 - 1 -. Note that by virtue of the self- 
consistency only phases commensurate with the lattice 
structure can be stabilized by a DMFT calculation, but 
not a truly incommensurate phase. Instead, we rather 
observe oscillations in the magnetization and occupation 
during the DMFT self-consistency cycle. Together with 
evidence from the one orbital model we can conclude that 
these oscillations, appearing when doping the Neel state 
at half filling, arc indeed signs of an incommensurate spin 
density wav o 33 i 34 ' 35 . For occupations of (n) < 1.5 the 
physical situation starts to become influenced by the dou- 
ble exchange mechanis m 36 i 37 and we thus expect to find 
ferromagnetic order- 3 ^^^^^^^ . Figure Q] and [2] 
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FIG. 2: (color online) The same as in figure[T]for J = 1.5VF. 
Compared to J = W the ferromagnetic phase extends to 
higher temperatures and larger hole doping. 

display the ferromagnetic polarization versus the occu- 
pation for two values of the Hund's coupling J for fixed 
U = 4W and U' = U — 2 J. One can see how the ferro- 
magnetic phase is stabilized at fillings n < 1.5. The blank 
rectangles close to half filling in figures [1] and [5] represent 
the already mentioned parameter regimes where we found 
no convergent solution to the DMFT, but an oscillatory 
behavior we interpreted as spin density waves. As they 
cannot really be stabilized, the true phase boundaries 
cannot be determined in this parameter region. 

As the ferromagnetic double exchange is mostly due to 
the Hund's coupling, it is not surprising that for increas- 
ing Hund's coupling the ferromagnetic state is stabilized 
up to higher temperatures and larger doping and can 
even extend beyond quarter filling. However, one must 
emphasize that Hund's coupling alone is not sufficient to 
enforce ferromagnetic order in the two orbital model. We 
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FIG. 3: (color online) Ferromagnetic order at quarter filling 
for T = and U = 4W as function of U' and J. Tuning 
Hund's coupling and the inter-orbital density-density interac- 
tion, one can observe a transition between an orbitally ordered 
ferromagnetic insulator and orbitally degenerate ferromag- 
netic metal. The symbols denote the parameters for which 
calculations were actually performed. The phase boundaries 
are fits to the calculated points and meant as guide to the 
eye. 



furthermore observe that one in addition needs a rather 
strong interaction parameter U to stabilize extended re- 
gions of ferromagnctism. For example, for U /W — 2 we 
found no ferromagnctism for 1 < n < 2. 



B. Magnetic phases at quarter filling, (n) = 1 

Quarter filling, like half filling, represents a very special 
point for the two orbital Hubbard model. In a classical 
picture there is one electron per site, which can choose 
between two orbitals. This picture makes already clear 
that orbital degeneracy and fluctuations will play an im- 
portant role at quarter filling. In Fig. [3] we show the 
ferromagnetic ground state phase diagram for n = 1 as 
function of U 1 and J. As noted before, for strong Hund's 
coupling J and moderate U' we obtain the orbitally ho- 
mogeneous ferromagnetic phase discussed in figure [U 
However, for large enough repulsive inter-orbital density- 
density interaction U' , we observe that the ferromagnetic 
spin alignment is accompanied by an antiferro-orbital or- 
der of the conventional Neel-typ o 39 i 45 i 46 i 47 . There is a 
first order transition between the orbitally ordered ferro- 
magnetic state and the homogeneous one: both the mag- 
netic and orbital polarization show jumps when crossing 
the phase boundary, see figured] Note that both states 
are strongly polarized, i.e. the jump in the magnetization 
is comparatively small. A more important aspect is that 
the orbitally ordered ferromagnetic phase is an insulator, 
while the homogeneous one is a metal. This particular 
difference becomes apparent from Fig.0 which shows the 
spectral function of both ferromagnetic states at quarter 
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FIG. 4: (color online) Orbital occupation (n m ,"r) ( m denotes 
the orbital) and spin polarization for U/W = 4, J/W = 1, 
T = 0, and different interaction strengths 17' . At {7'/W ~ 1.4 
the transition between the homogeneous ferromagnetic state 
(f/'/W 7 < 1-4) to the orbital ordered ferromagnet (U' /W > 
1.4) occurs. Lines are meant as guide to the eye. 



filling. The left (right) panel illustrates the spin and or- 
bital configuration of the metallic (insulating) phase in 
the upper part and the spectral function corresponding 
to that phase in the lower part. The Fermi energy lies 
at uj = 0. As can be easily seen, the homogeneous ferro- 
magnetic state has a large spectral weight at the Fermi 
energy, thus representing a metallic system. On the other 
hand, the orbitally ordered ferromagnetic state has a gap 
around the Fermi energy, thus representing an insulating 
state. The gap width of the insulating state decreases 
when approaching the transition line. Thus, by vary- 
ing Hund's coupling J or the inter-orbital interaction U' 
one can observe a metal insulator transition (MIT) be- 
tween two almost fully polarized ferromagnetic phases. 
Note that this metal insulator transition is very different 
from the usual paramagnetic one, which appears in the 
Hubbard model at half filling as function of U. The MIT 
observed here between the ferromagnetic phases is rather 
due to a strong inter-orbital density-density interaction, 
which is responsible for driving the orbital ordering. As 
the latter introduces a doubling of the unit cell, the in- 
sulating solution is thus akin to the antifcrromagnctic 
insulator. The usual intra-orbital Hubbard interaction 
U plays a minor role in this transition. 

It is worth noting that within a simplified two-site 
model representing the AB structure of the Neel-state 
the antiferro-orbital situation wins over the orbitally ho- 
mogeneous for all values of U' and J due to virtual hop- 
ping and the gain of Hund's exchange energy. It is thus 
the presence of the lattice which allows for smaller U' 
through an additional gain in kinetic energy the forma- 
tion of the homogeneous metallic ferromagnet. In view of 
effects like colossal magneto resistivity it is an obviously 
interesting question, what influence external control pa- 
rameters like temperature and magnetic field will have 
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FIG. 5: (color online) Spectral function and sketch of the elec- 
tronic configuration. Left panel: metallic ferromagnetic phase 
for U/W = 4, J/W = 1.5 and U' = U - 23. Right panel: 
insulating ferromagnetic phase for U /W = 4 and J/W = 0.5. 
uj — represents the Fermi energy. 
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FIG. 6: Orbital ordered ferromagnetic insulator for U = 4W, 
J — W/2, U = U — 2 J as function of the chemical potential. 
The upper (lower) panel shows the occupation (polarization). 
Notice the jump in these quantities. 



when one is close to the phase transition. These ques- 
tions are presently under investigation. 

Another important aspect is the dependence of the 
phases upon doping. For the homogeneous ferromag- 
netic phase we already saw when discussing figure [21 that 
its filling can be varied smoothly. The dependency on 
the chemical potential of the orbital ordered ferromag- 
net, on the other hand, is dramatically different and can 
be seen in figure El The plot shows that for a critical 
chemical potential the filling of the system jumps from 
approximately n = 0.88 to quarter filling n = 1. At 



5 



the same critical chemical potential also the polarization 
jumps from a non-polarized phase to nearly fully polar- 
ized. In other words, there is a first order transition 
between a paramagnetic phase with filling less than one 
and an orbitally ordered ferromagnetic state at quarter 
filling. Consequently the electronic system shows phase 
separation between these states and the precise physics 
will depend on additional interactions, like long-range 
Coulomb interaction, or additional degrees of freedom 
like the lattice. As already noted before, this appear- 
ance of phase separation is similar to what we find in 
the Hubbard model in the antiferromagnetic phase at 
half fillin g 31 i 48 ' 49 i 50 , and thus seems to be a generic fea- 
ture of the symmetry broken phases on an AB lattice. 
Note, however, that longer-ranged hopping can actually 
destroy this phase separation^, depending on the sign of 
the additional hopping and the type of doping. 

In the region of large U' , where the orbitally ordered 
ferromagnet is found, we were also able to stabilize an 
antiferromagnetic phase at quarter filling. As the sys- 
tem at quarter filling is dominated by ferromagnetic dou- 
ble exchange, one actually does not expect such a phase 
here. Like the orbitally ordered ferromagnet, this anti- 
ferromagnetic phase also exists only exactly at quarter 
filling. The spectral function and the doping dependence 
can be seen in figure [7J The spectral function shows 
that also this state is a perfect insulator. When trying 
to dope it we again find phase separation to a paramag- 
netic metal away from quarter filling. In order to find out 
which state is the thermodynamically stable one, we cal- 
culated the energy of both states. The result is that, as 
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FIG. 7: (color online) Antiferromagnetic insulator for U = 
4W, J = W/2, U' = U — 2Jas function of the chemical po- 
tential. The upper (lower) panel shows the occupation (polar- 
ization). Notice the jump in the properties. The right panel 
shows the spectral function. 

expected, the orbital ferromagnet has the lower energy, 
thus is the thermodynamically stable one. Looking at 
the different terms in the energy, the kinetic energy gives 
a larger decrease for the antiferromagnetic state than for 
the ferromagnetic state, while the interaction terms in- 



crease the energy of the antiferromagnet. Varying the 
parameters J and U', we always find the antiferromag- 
netic state having the higher energy. 

C. Charge ordering at quarter Filling 

For large Hund's coupling J/2 > U', the term U' — J/2 
defining the inter-orbital density-density interaction be- 
comes attractive. Although at first glance such a large 
Hund's coupling appears unphysical, one might have sit- 
uations, for example Jahn- Teller coupling to phonons, 
which can lead to additional contributions to U' , typ- 
ically reducing it effectively In this case such an at- 
tractive interaction can effectively be generated and the 
physics will change dramatically. The first thing one 
notes is that during a numerical calculation it becomes 
very difficult to stabilize fillings other than n = 2 or 
n = 0. Inspired by this difficulty we investigated charge 
ordered phases in this parameter regime; and for suffi- 
ciently large Hund's coupling J it is indeed possible to 
stabilize a charge ordered state with alternating almost 
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FIG. 8: Occupation of neighboring sites in the charge ordered 
state for J/W = 2. The inter-orbital interaction becomes 
repulsive for U' /W > 1. The lines are meant as guide to the 
eye. 

doubly occupied and nearly empty sites and an average 
occupation of n = 1, sec figure [H Interestingly, this 
state seems to reach into the regime with U' > J/2, al- 
ready representing a repulsive inter-orbital interaction. 
If V is increased further, this charge ordered state fi- 
nally becomes unstable for U'JW ~ 1.7, i.e. we obtain a 
quarter filled paramagnetic state. Note that we do not 
observe a vanishing order parameter but rather a jump 
as U' y U' c , indicating a first order transition here, too. 
Another open point is how the charge ordered state con- 
nects to the magnetic phases present in this parameter 
region. As we have not yet been able to perform calcu- 
lations allowing for both charge and magnetic order, we 
cannot tell whether there is a direct transition into one 
of the ferromagnetic phases rather than to the param- 
agnetic state. For the magnetic properties of the charge 
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ordered state one may expect some kind of magnetic or- 
der between the half filled sites, which however requires 
larger unit cells to be used in the calculations. Besides 
magnetism, the charge density wave state is a perfect in- 
sulator, too, and again it is not possible to dope this state 
away from quarter filling. 

IV. SUMMARY 

In conclusion, we have analyzed the magnetic phase 
diagram of a two orbital Hubbard model within DMFT 
and NRG. While around half filling the system behaves 
quite similar to the one orbital Hubbard model, there 
occurs an extended ferromagnetic phase for occupation 
(n) < 1.5. Such a ferromagnetic phase cannot be ob- 
served for the one orbital case on a Bethe lattice with 
semi elliptic DOS and is due to double exchange mech- 
anism present in the two-orbital model. With increas- 
ing Hund's coupling J this ferromagnetic phase becomes 
more and more extended and finally can also be observed 
for occupation smaller than 1. 

A particularly important point in the phase diagram 
is quarter filling, where we could observe four different 
ordered phases, and an especially interesting feature of 
the quarter-filled case is the presence of a metal-insulator 
transition between two ferromagnetic states. The transi- 
tion seems to be of first order and is driven by the inter- 
orbital density-density interaction U' . For low U' the 
ferromagnetic state is homogeneous and metallic. For 
large U' , an orbital order can be observed in addition to 



the ferromagnetic one, which now is accompanied by an 
insulating behavior as the orbital order breaks transla- 
tional symmetry similar to the Nccl state at half filling. 

Besides these ferromagnetic states, we also could sta- 
bilize an antiferromagnetic insulating and a charge or- 
dered insulating state at quarter filling. The antiferro- 
magnetic state exists in the same parameter region as 
the orbital ordered ferromagnetic state, but has a higher 
energy, i.e. will not be the thcrmodynamically stable one. 
The charge ordered state can be observed for rather large 
Hund's coupling J respectively reduced inter-orbital den- 
sity interaction U' . The latter situation can for example 
be realized in the presence of Jahn- Teller phonons, which 
definitely play an important role in two-orbital systems 
with e g symmetry. Thus, the inclusion of lattice degrees 
of freedom is a very important extension and presently 
under investigation. 
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